A hyperbolic slicing condition adapted to Killing fields and densitized lapses 
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We study the properties of a modified version of the Bona-Masso family of hyperbolic slicing 
conditions. This modified slicing condition has two very important features: In the first place, it 
guarantees that if a spacetime is static or stationary, and one starts the evolution in a coordinate 
system in which the metric coefficients are already time independent, then they will remain time in- 
dependent during the subsequent evolution, i. e. the lapse will not evolve and will therefore not drive 
the time lines away from the Killing direction. Second, the modified condition is naturally adapted 
to the use of a densitized lapse as a fundamental variable, which in turn makes it a good candi- 
date for a dynamic slicing condition that can be used in conjunction with some recently proposed 
hyperbolic reformulations of the Einstein evolution equations. 
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I. INTRODUCTION 



Specifying a good foliation of spacetime is of funda- 
mental importance when studying the dynamical evo- 
lution of systems with strong gravitational fields. In 
the 3+1 decomposition of General Relativity, a foliation 
of spacetime into spacelike hypersurfaces is described in 
terms of the lapse function a that measures the interval 
of proper time between neighboring hypersurfaces along 
their normal direction. The choice of a particular folia- 
tion represents the freedom one has to specify the time 
coordinate and is therefore arbitrary. In practice, how- 
ever, one can not choose a foliation ahead of time since 
one is trying to solve for the geometry of the spacetime 
itself, so one must choose instead a "slicing condition" , 
that is, some geometric or algebraic condition that allows 
one to calculate the lapse function dynamically during 
the evolution. Slicing conditions come in many different 
forms, and they are usually chosen by balancing their 
ease of implementation with the need to obtain a well 
behaved coordinate system, as well as the well-posedness 
of the resulting system of evolution equations. 

In a recent paper [1], one of us studied a particular 
slicing condition that has the property of being obtained 
through a hyperbolic evolution equation for the lapse, 
and in many cases allows one to construct a foliation that 
avoids different types of pathological behaviors. This slic- 
ing condition, known as the Bona-Masso (BM) slicing 
condition [2], has been used successfully for many nu- 
merical simulations of strongly gravitating systems such 
as black holes (see for example [3, 4, 5]). Other forms 
of hyperbolic slicing conditions, which generalize in one 
way or another the original proposal of Bona and Masso, 
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have been suggested in the last few years, examples of 
which can be found in [5, 6], and most recently in [7]. 



In this paper we argue that the original form of the 
BM condition (as well as most of its generalization) has 
two important drawbacks: First, it is not well adapted 
to the presence of Killing fields in the sense that if one 
uses it to evolve a static spacetime, it can easily drive the 
time lines away from the Killing direction. Second, the 
standard BM condition is not well adapted to the case 
when one wants to use a densitized lapse as a fundamen- 
tal variable. Such a densitized lapse has been recently 
advocated in the context of hyperbolic reformulations of 
the Einstein equations (see for example [8, 9]), and is 
therefore an important issue to consider. Here we will 
study a modified version of the BM slicing condition that 
addresses both these issues at the same time. This modi- 
fied BM condition has already been used in the literature, 
but as far as we are aware its properties have never been 
studied in any detail. On this manuscript we will limit 
ourselves to studying the properties of this slicing condi- 
tion independently of the Einstein equations, and leave 
the study of how this condition couples to the Einstein 
evolution equations for a separate work [10]. 



This paper is organized as follows. In Sec. II we make a 
brief introduction to the BM family of slicing conditions. 
Section III motivates the introduction of a modified BM 
condition from the point of view of compatibility with 
static and stationary solutions. We then study the rela- 
tion of the modified condition both with densitized lapses 
and with the divergence of the coordinate time lines, and 
introduce a coordinate independent way of writing the 
condition. In Sec. IV we analyze the hyperbolicity of the 
modified BM condition. We start by discussing the way 
in which the shift vector is chosen, and later analyze the 
hyperbolicity of the coupled slicing-shift evolution sys- 
tem. Finally, Sections V and VI study under which cir- 
cumstances the modified BM condition avoids focusing 
singularities and gauge shocks. We conclude in Sec. VII. 
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II. THE BONA-MASSO FAMILY OF SLICING 
CONDITIONS 

The BM family of slicing conditions [2] has been dis- 
cussed extensively in the literature (for a detailed dis- 
cussion see [11] and references therein). Here we will 
limit ourselves to making a very brief description of its 
main properties. The BM slicing condition is obtained 
by asking for the lapse function to satisfy the following 
evolution equation 

j t a=(d t -jCfs)a = -a 2 f(a)K, (2.1) 

with C/3 the Lie derivative with respect to the shift vec- 
tor /3 l , K the trace of the extrinsic curvature and f(a) a 
positive but otherwise arbitrary function of a. The con- 
dition above is a generalization of slicing conditions that 
have been used in evolution codes based on the Arnowitt- 
Deser-Misner (ADM) formulation [12, 13] since the early 
90's [14, 15], and was originally proposed in the context 
of the Bona-Masso hyperbolic re-formulation of the Ein- 
stein equations [2, 16, 17, 18, 19]. It is, however, very 
general and can be used with any form of the 3+1 evolu- 
tion equations. Since we will later introduce a modified 
version of the BM slicing condition, we will refer from 
now on to condition (2.1) as the "standard BM condi- 
tion" . 

A very important property of the standard BM con- 
dition is the fact that the shift terms included through 
the Lie derivative in equation (2.1) are such that one is 
guaranteed to obtain precisely the same spacetime folia- 
tion regardless of the value of the shift vector. At first 
sight, this would seem to be a natural requirement for 
any slicing condition, but as we will see below this might 
not be the most important property a slicing condition 
must have. 

The BM slicing condition can be shown to lead to a 
generalized wave equation for the lapse 



-a 3 f [KijK* - (2/ + af) K 2 ] . (2.2) 

From this equation one can easily see that the wave speed 
along a specific direction x l given by 

v g =ay/W\ (2.3) 

where jij is the spatial metric. The above expression ex- 
plains the need for f(a) to be positive: If it weren't, the 
wave speed would not be real and the equation would be 
elliptic instead of hyperbolic. Notice also that the gauge 
speed (2.3) can be smaller or larger that the speed of light 
depending on the value of /. Contrary to what one might 
expect, having a gauge speed that is larger than the speed 
of light does not lead to any causality violations, as the 
supcrluminal speed is only related with the propagation 
of the coordinate system. In fact, empirically the most 



successful slicing conditions for the simulation of black 
hole spacetimes have been precisely those that allow su- 
pcrluminal gauge speeds (an example of this is the 1+log 
slicing condition mentioned below). 

Reference [1] shows that, for some specific choices 
of the function f(a), the standard BM slicing condi- 
tion can avoid both focusing singularities [19] and gauge 
shocks [20]. We will expand on these points in the fol- 
lowing sections. 

A very important particular case of the BM condition 
corresponds to the choice /(a) = 1, which leads to the so- 
called "harmonic slicing" , for which the time coordinate 
obeys the simple wave equation: 

□* = 9^K» = . ( 2 - 4 ) 

with g^ v the spacetime metric. That is, / = 1 corre- 
sponds to the case when the time coordinate is a har- 
monic function. One can easily show that for harmonic 
slicing one also has the following relation 

a = htf) 7i /2 , (2.5) 

with h{x l ) an arbitrary time independent function and 
V / 7^T the volume element associated with observers mov- 
ing normal to the hypcrsurfaces (which will differ from 
the coordinate volume element in the case of a non- 
vanishing shift vector). 

Another particular case worth mentioning corresponds 
to "1+log" family, for which f(a) = N/a. In this case 
we find the following relation between the lapse function 
and the normal volume elements 

a = h(x i )+ln{ 1 ^ 2 ) , (2.6) 

which explains the name 1+log. The case N — 2 has been 
found empirically to be particularly well behaved [3, 4, 5], 
something which Ref. [1] attributes (a posteriori) to the 
fact that this is the only member of the 1+log family that 
avoids gauge shocks even approximately. Notice that for 
N = 2 the gauge speed in asymptotically flat regions 
where a ~ 1 becomes \[2 > 1. In regions inside a black 
hole, where the lapse typically collapses to zero, the gauge 
speed can become extremely large. 

Having briefly discussed the main properties of the 
standard BM slicing condition, we will now turn to our 
proposal for a modified version of this condition, moti- 
vated by the analysis of stationary spacetimes. 

III. A MODIFIED BONA-MASSO SLICING 
CONDITION 

Let us consider for a moment a static or stationary 
spacetime, and let us assume that we have chosen a co- 
ordinate system in which the metric coefficients are al- 
ready time independent. Further, we will assume that 
our coordinate system is such that the shift vector does 
not vanish. One would think at first sight that having 
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a non- vanishing shift implies that the spacetime can not 
be truly static but is at most only stationary. This is not 
entirely correct, as one can write a static spacetime in 
a coordinate system with non-vanishing shift and where 
the metric coefficients are nevertheless still time inde- 
pendent. An example of this is the Schwarzschild metric 
written in Kerr-Schild or Painleve-Gullstand coordinates. 

Assuming we have such a situation, the ADM evolution 
equations imply in particular that 



dtHj = -2aKij + Cpjij 

= -2aKij + Dipj + DjPi = , 



(3.1) 



with Di the spatial covariant derivative. From the last 
equation one can obtain 



(3.2) 



with 7 1 / 2 the coordinate volume elements. The trace of 
the extrinsic curvature must therefore be given by 



K = 



(3.3) 



Let us now see if the standard Bona-Masso slicing con- 
dition is compatible with the time-independent charac- 
ter of our spacetime in the sense of predicting a time- 
independent lapse as well. Substituting our expression 
for K into equation (2.1) we find 



d t a — [3 l diOt 



2a 2 f(a)K 
2af(a)D i f3 i 



(3.4) 
(3.5) 



It is clear from looking at this expression that for arbi- 
trary /(a), dta will generally not vanish. That is, the 
lapse will evolve away from its preferred value and as a 
consequence the spatial metric will not remain time in- 
dependent. This indicates that the standard BM slicing 
condition is not well adapted to the evolution of station- 
ary spacetimes. 

A modified version of the BM slicing condition that is 
well adapted to such time independent spacetimes can 
nevertheless be easily obtained by just asking for 



d t a 



af(a) 



7 



1/2 



da 



1/2 



(3.6) 



which automatically guarantees that the lapse will not 
evolve if the spatial metric is time independent. One can 
integrate the last equation trivially to find 
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x ' 2 = F{x l ) exp 



da 
af{a) 



(3.7) 



with F{x l ) a time independent function. This tells us 
that there is a very general functional relationship be- 
tween a and 7 1 ' 2 . For f(a) = constant this relation- 
ship is a power law, but in other cases it is more general 
(for example, the well known "1+log" slicing is usually 
obtained from / = 2/ a, which gives us an exponential 



relation between a and 7 1 / 2 ). Functional relationships 
between a and 7 1 / 2 have been considered before in the 
context of finding hyperbolic re-formulations of the 3+1 
evolution equations. For example, in reference [21], Frit- 
telli and Reula propose a general power law relation, 
which as just mentioned is a particular form of (3.7) 
with f(a) constant. More recently, Sarbach and Tiglio [6] 
have considered a completely general functional relation- 
ship with the sole restriction that da/d-f > 0, which 
clearly includes (3.7) (asking for f(a) > guarantees 
that da/dj > 0). What makes (3.7) more interesting 
than a very general relationship is the fact that one can 
learn a lot about the properties of the slicing by studying 
the effect of different forms of f(a). 

Substituting now (3.2) into (3.6) we find 



d t a = -af(a) (aK - D^) 



(3.8) 



This is the modified version of the BM slicing condition 
we wish to study here (compare with (2.1)). Several com- 
ments are in order here. First, we should mention that 
the modified BM condition defined above has been used 
before in the literature in Refs. [22, 23]. Both these ref- 
erences study the numerical evolution of a Schwarzschild 
spacetime written in Kerr-Schild coordinates, and at- 
tempt to maintain the static (exact) solution stable dur- 
ing the numerical simulation. Because of this they need 
to use a slicing condition that maintains the lapse equal 
to its initial value in the continuum limit, but allows some 
dynamics to respond to numerical truncation errors at fi- 
nite resolutions. The modified BM condition is used in 
those references precisely for the reason outlined above, 
but no suggestion is made to use it in the more general 
case or to study its properties. 

Second, as condition (3.8) does not include the Lie 
derivative terms of the lapse with respect to the shift vec- 
tor, it will not give us the same foliation of spacetime for 
a different choice of the shift vector, i.e. the foliation of 
spacetime one obtains will depend on the choice of shift. 
We believe this is not a serious drawback since in a par- 
ticular situation one would presumably want to choose a 
slicing condition and a shift vector that are closely inter- 
related. 

Also, we believe that using a slicing condition that is 
compatible with a static solution is a necessary require- 
ment if one looks for symmetry seeking coordinates of 
the type discussed by Gundlach and Garfinkle [24] and 
by Brady et.al [25], that will be able to find the Killing 
fields that static (or stationary) spacetimes have, or the 
approximate Killing fields that many interesting astro- 
physical systems will have at late times. Of course, hav- 
ing a slicing condition that is compatible with staticity is 
not enough, one also needs to have a shift condition that 
has the same property. Otherwise, the shift evolution will 
also drive us away from the frame in which the staticity 
is apparent. In this paper we will not deal with the issue 
of the shift choice, but will consider it in a future work. 

And finally, even if one has gauge conditions (for both 
lapse and shift) that keep a static solution static, this 



4 



does not guarantee that one can find the Killing direction 
if one starts in the wrong coordinate system. All that it 
guarantees is that if we do find this Killing direction we 
won't be driven away from it again. 



A. Densitized lapse 

Since the early 90's, many re-formulations of the 3+1 
evolution equations of general relativity have been pro- 
posed [2, 8, 9, 17, 18, 19, 21, 26, 27, 28, 29, 30, 31, 32, 33]. 
The main purpose of most such re-formulations has been 
to recast the Einstein evolution equations as a strongly or 
symmetric hyperbolic system motivated by the fact that 
for such systems one can prove that the evolution equa- 
tions are well posed, and with the hope that such well- 
posedness will bring with it improvements in numerical 
simulations, both in terms of stability and in terms of al- 
lowing one to impose boundary conditions in a consistent 
and well behaved way. 

While it is not the purpose of this paper to study the 
different hyperbolic formulations of the Einstein evolu- 
tion equations, there is one related point that is crucial 
for the discussion of slicing conditions. It was realized 
early on that there were two main problems when try- 
ing to go from the standard ADM evolution equations 
to a strongly or symmetric hyperbolic formulation. One 
problem was associated with the existence of constraint 
violating solutions that spoil hyperbolicity and required 
for its solution the addition of multiples of the constraint 
equations, in a number of different ways, to the evolution 
equations (the constraint violating modes are not elimi- 
nated, just transformed in a way that allows a hyperbolic 
system to be constructed). 

The second problem, more directly related to our dis- 
cussion, was the observation that it is not possible to 
construct a strongly hyperbolic formulation of the 3+1 
evolution equations if the lapse function is considered to 
be an a priori known function of space and time. Two 
different routes have been followed to solve this second 
problem. The first route, taken in the Bona-Masso hy- 
perbolic re-formulation [2, 17, 18, 19], was to propose an 
evolution equation for the lapse (equation (2.1) of the 
previous section), and then construct a strongly hyper- 
bolic system of equations where the lapse was considered 
just another dynamical variable. This same route has 
been followed very recently by Lindblom and Scheel [7], 
where generalizations of standard BM slicing condition 
and the 'T-driver" shift condition [5] have been used to 
construct a symmetric hyperbolic system that includes 
lapse and shift as dynamical variables. The second route 
has been to take not the lapse, but rather the densitized 
lapse q := cr^ 1 / 2 , to be a prescribed function of space 
and time (see [8, 9] and references therein). Both these 
routes have been successful in constructing strongly hy- 
perbolic re-formulations of the Einstein evolution equa- 
tions. The different approaches are related, but are not 
equivalent, as can be seen easily if we consider for a mo- 



ment the BM slicing condition in the case / = 1 (the 
harmonic slicing case). As we have seen, in that case the 
lapse takes the form 



a 



(3.9) 



1 /2 

with 7„ the volume elements associated with the nor- 
mal observers, and h(x l ) an arbitrary function of space. 
We can in fact turn the function h(x % ) into an arbitrary 
function of both space and time if we add a source term 
to the standard BM condition in the following way 



{d t -Cp)a= -a 2 f(a) K + H(x, t) 



(3.10) 



with H (x, t) arbitrary. On the other hand, the condi- 
tion for the densitized lapse to be a known function of 
spacetime takes the form 



a = q(x,t) 7 



1/2 



(3.11) 



where 7 1 / 2 are the coordinate volume elements. In the 
case of a vanishing shift vector, normal and coordinate 
volume elements coincide, and the BM condition can be 
seen as a generalization of the prescribed densitized lapse 
condition. But for non-zero shift vector this is no longer 
the case. 

One can in fact rewrite the standard BM condition in 
terms of the densitized lapse q in the following way 



(dt-Cp) q = -q 2 1 1 ' 2 (f-l)K, 



(3.12) 



which shows that even for / = 1, the densitized lapse q 
will evolve dynamically driven by the Lie derivative term. 

The crucial observation here is that the modified BM 
condition (3.8), when written in terms of the densitized 
lapse, takes the form 



d t q = -q(f-l) \qi 1/2 K-D t p 



(3.13) 



It is now clear that by taking / = 1 this equation re- 
duces to the case of a static densitized lapse (the case of 
a prescribed densitized lapse that is not time indepen- 
dent can be easily obtained by adding a source term to 
the above equation). The modified BM slicing condition 
would therefore seem to be a natural generalization of the 
prescribed densitized lapse, and should be well adapted 
to hyperbolic re-formulations of the Einstein equations 
that use of a densitized lapse [10]. 



B. Divergence of the time lines and Killing fields 

From the ADM equations one can easily show that 
the divergence of the time lines is given by the following 
expression 



V^f = - [d t a - a (aK - Dtf 1 



(3.14) 
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where is the vector tangent to the 3+1 time lines, 
which is defined in terms of the normal and shift vectors 
as 

V = an 11 + . (3.15) 

Equation (3.14) implies that the modified BM slicing 
condition can also be written as 

d t a = a V„f • (3-16) 

The last equation shows that the evolution of the lapse 
is directly related to the divergence of the time lines. We 
then see that just as the original BM slicing condition 
was such that the lapse reacted to the divergence of the 
normal observers, the modified condition ensures that 
the lapse reacts to the divergence of the coordinate time 
lines. 

Let us assume for the moment that our spacetime has 
a future-pointing Killing field v+ In that case we will 
have 

V„w v + V v u M = V/ = 0. (3.17) 

If we now assume that our time lines are oriented along 
the Killing direction, then equation (3.16) automatically 
implies that the lapse function will be time independent. 
This, of course, we already knew since it was our initial 
motivation for modifying the BM condition. 

C. Generalized wave equation for the time function 

In reference [1] it was shown that the standard BM 
slicing condition can be written as a generalized wave 
equation for a "time function" <p m the following way 

(r-anV)V^ = 0, (3.18) 

with the unit normal vector to the spatial hypersur- 
faces and a := 1/ f(a) — 1. The different members of the 
spacetime foliation can then be obtained as the level sets 
of the time function <\>. 

One can also construct such a foliation equation for the 
modified BM condition (3.8). The corresponding equa- 
tion for the time function <j> turns out to be 

(g^ - V nj V M V„ = V M (^j . (3.19) 

By writing equation (3.19) in the standard 3+1 coor- 
dinate system adapted to the foliation, it is not difficult 
to show that it is in fact equivalent to the slicing condi- 
tion (3.8). On the other hand, notice that when written 
in a different coordinate system, the vector does not 
have to be aligned with the time lines any longer. Also, 
the shift vector /3 M will have a non-zero time component. 
However, we will still have 

/3"n„ = . (3.20) 



In the case of vanishing shift, equation (3.19) reduces 
to equation (3.18), which is just another way of saying 
that in that case the original and modified BM conditions 
coincide. For non- vanishing shift, however, both equa- 
tions differ. The foliation equation for the original BM 
condition makes no reference to the shift, which implies 
that the foliation is independent of our choice of shift (as 
already mentioned above). The modified foliation equa- 
tion, however, clearly depends on the shift choice, so the 
foliation of spacetime one obtains will depend on the shift 
as well. 

The foliation equation (3.19) is very useful when trying 
to understand the properties of the slicing condition in 
a covariant way that is independent of the Einstein field 
equations. 



IV. HYPERBOLICITY 

The concept of hyperbolicity is of fundamental impor- 
tance in the study of the evolution equations associated 
with a Cauchy problem. Some measure of hyperbolicity, 
even in a weak sense, implies that the system of equa- 
tions is causal, i.e. that the solution at a given point 
in spacetime depends only on data in a region of com- 
pact support to the past of that point (the characteristic 
cone). Stronger versions of hyperbolicity can also be used 
to prove rigorously that the system of equations is well- 
posed, that is, that its solutions exist (at least locally), 
are unique, and are stable in the sense that small changes 
in the initial data will correspond to small changes in the 
solution. Hyperbolicity also allows one to construct well- 
posed initial-boundary problems, which implies that one 
should be able to obtain well behaved boundary condi- 
tions for numerical simulations with artificial boundaries. 

Because of this, showing that a given system of evo- 
lution equations is hyperbolic has become an important 
test for new formulations of the 3+1 equations. In our 
case, since we are studying a slicing condition that is ob- 
tained through an evolution equation for the lapse, the 
question of the hyperbolicity of the gauge condition be- 
comes important. Since we want to look at this issue in 
a way that is independent of the Einstein equations, we 
will consider from now on a given background spacetime 
(which may or may not obey Einstein's equations), and 
study our slicing condition on this fixed background. 



A. Prescribed shift vector 

Since the foliation equation (3.19) involves the shift 
vector /?'*, in order to analyze its hyperbolicity we must 
say something about how the shift vector evolves in time. 
The simplest approach is to assume that we have a pre- 
scribed, i.e. non-dynamical, shift vector. However, this 
immediately leads us into the question: What does it 
mean to have a prescribed shift vector in a covariant 
sense? 
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Clearly, we can not just assume that the shift vector 
is an a priori known vector-field in spacetime. This is 
because the shift vector must always be parallel to the 
spatial hypersurfaces making up our foliation, and those 
hypersurfaces are precisely what we are trying to solve 
for. This means that inevitably, as we solve for the hy- 
persurfaces, the shift vector must evolve dynamically to 
guarantee that it remains parallel to them. At most, we 
can ask for the magnitude and direction of the shift vec- 
tor within a given hypersurface to be prescribed functions 
of space and time. 

We will then propose the most general evolution equa- 
tion for the shift that is compatible with the fact that the 
shift lives on the spatial hypersurfaces. For this we pro- 
ceed as follows: On a given hypersurface we can choose 
a basis of spatial vectors ef , and express the shift vector 
in terms of such a basis: 



/3" = ef b* 



(4.1) 



where the b 1 are the components of the shift vector in the 
basis under consideration. We will now identify this ba- 
sis with the standard 3+1 spatial coordinate basis. This 
means that the spatial basis ef , together with the time 
vector form a coordinate basis for the spacetime at 
the point under study (they form the standard 3+1 coor- 
dinate basis), which implies the following commutation 
relation 



f d v e? = e» d v V . 



(4.2) 



Using this relation, we can obtain the following very 
general evolution equation for the shift vector: 



(4.3) 



where s M = 1" d v V . Notice that the last equation is 
fully covariant even if it is written in terms of partial 
derivatives, as the Christoffel symbols cancel out. If we 
now assume that the components V of the shift in the 
3+1 spatial coordinate basis are prescribed functions of 
spacetime, then we can consider the s M as source terms 
in the hyperbolicity analysis. 

We will use equation (4.3) above as our general evolu- 
tion equation for the shift. Notice that in the particular 
case when we restrict ourselves to the 3+1 coordinate 
system, for which f = (1,0,0,0), this equation simply 
reduces to 

dtp 1 =d t b\ (4.4) 
which is of course to be expected. 

B. Hyperbolicity of the foliation equation 

We are interested in studying the hyperbolicity of the 
modified BM slicing condition in a way that is indepen- 
dent of the Einstein equation. In order to do this we 



assume we have a fixed background spacetime with met- 
ric gfiv, and study the foliation equation (3.19), which we 
repeat here for clarity: 



(4.5) 



Here ^ is a scalar function whose level sets identify the 
elements of the foliation, and a := 1/ f(a) — 1. In terms 
of 4>, the lapse function and the unit normal vector can 
be expressed as 



a = (-V M 0V. 



-1/2 



n M = -aV M . 



(4.6) 
(4.7) 



We will now concentrate on a point on a given slice, and 
construct locally flat coordinates in its neighborhood. We 
will further define the following first order quantities 



dt<f> = n , 

di<j) = . 

The lapse and unit normal vector then become 

T 2 ,t,2\-V2 



a = (n 2 -* 2 y 

= a(II,-* J ) 



(4.8) 
(4.9) 



(4.10) 
(4.11) 



with "J 2 = ^2 "f 2 . Notice that, as we are in locally flat co- 
ordinates, there is no difference between lower and upper 
spatial indexes, so we will be using them indiscriminately. 

On the other hand, since the shift vector is parallel to 
the hypersurface we must have 

/?"n M = , (4.12) 

which allows us to express the /3° component as 







n 



(4.13) 



With the definitions above we can now rewrite the fo- 
liation equation (4.5) as: 



- (Pa 2 - ((1 + a)n 2 + * 2 ) Q) d t u 



a 2 n 



dt(3 l 



+B m d rn n + c mj d m ^j - — d m p m = o , (4.14) 



where dot stands for the partial time derivative, and 
where we have defined 



P := (l + a)n 2 -* 2 



Q := 



IP ' 



(4.15) 
(4.16) 



am 

B m := -((l + a)n 2 + * 2 ) tL- 



+ (2aa 2 - (2 + a)Q)n* m , (4.17) 
C mj ._ s m i -V m (aa 2 W -{2 + a)p j ) . (4.18) 
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Notice that to arrive at Eq. (4.14) above we have used 
the fact that 



(4.19) 



We can also rewrite the evolution equation for the shift 
introduced in the last section, Eq. (4.3), in our locally flat 
coordinates to find 

2a 2 * m /3 m * 4 a 4 n + iid t p l - N m d m u 

+L^ m d m ^j - V m d m l3 l - £l = , (4.20) 

where now 

N lm := nQ5 irn + 2a 2 n* 4 (Q* m + (i m ) , (4.21) 

L ijm ._ ( S ij + 2a 2 W&) I3 m . (4.22) 

Equations (4.14) and (4.20) can now be used to find 
the following evolution equations for II, and the spatial 
part of the shift vector fi l : 

o t u = i- (r 2 m 9 m n 

+ T£ m d m *j + T{ m d m P + F°) , (4.23) 
8t/? = Drd m n + D^ m d m ^ 3 



where 

Ti = 
T m = 



P(Q~a 2 ), 

-6 jm + a¥ (a 2 * m -/? m )- 



(4.24) 

(4.25) 
(4.26) 



/3 j * r 



i-.ijm 
u 2 



p.ijm 




1 / T-* m 

_ I 5H^m _ 2a 2 n 2g^i^_ 



(4.32) 
(4.33) 



Equations (4.23) and (4.24), together with Eq. (4.19), 
form our closed set of evolution equations. Notice that 
these equations are only valid if a 2 ^ Q, the case a 2 = Q 
being degenerate. The reason for this can be traced back 
to the fact that if a 2 = Q, then (3 zi f>i — (all) 2 , which in 



means that the vector i M has no time component and 
equation (4.3) stops being an evolution equation for the 
shift. That this is a purely coordinate problem can be 
seen from the fact that we can always boost our locally 
flat coordinates in such a way that becomes zero and 
the problem disappears. 

In order to determine whether or not the system of 
evolution equations is hyperbolic, we first write it in ma- 
trix notation. We start by defining the vector u in the 
following way 

u= (II, * X ,F, (3 y ,(3 z ) . (4.34) 

Thus the system of equations can be expressed as 

d t u = M x d x u + M y d y u + M z d z u + s , (4.35) 

where the Jacobian matrices M l and the source vector s 
depend on the w's but not their derivatives. 
The matrix 1VP has the particular form 



(4.36) 



with i, j — x, y, z. The matrices M y and M z have similar 
structures. 

Having written our system of equations in the 
form (4.35), we can now proceed to study its hyperbol- 
icity properties. In order to do this one should first con- 
struct the "principal symbol" M l tii, where rii is an arbi- 
trary unit vector. The system is then said to be strongly 
hyperbolic if the eigenvalues of the principal symbol are 
real and there is a complete set of eigenvectors for all n,; 
furthermore, the system is said to be symmetric hyper- 
bolic if the principal symbol can be symmetrized in a way 
that is independent of the rij [34]. If, on the other hand, 
the eigenvalues are real, but there is no complete set of 
eigenvectors, the system is only weakly hyperbolic. Here 
we will use a shortcut useful when all directions have 
equivalent structures: We will find the eigenvalues of one 
of the Jacobian matrices M l which by inspection can 
later be generalized to any arbitrary direction. We will 
concentrate on the matrix M x , as results for the other 
two matrices can be obtained afterward in a straightfor- 
ward way. 

The eigenvalues of M. x can be found to be 





/ T£/T1 


Ti x /Tl 


T{ x fT\ \ 




1 








M x = 
























\ Df 




DT ) 



A± 

A (3) 

A (4) 
A (6) 



^ (a* x ± R) , 



Il(Q-a 2 ) ' 

n ' 

A< 7 ) = , 



= A< 5 ) 



(4.37) 
(4.38) 

(4.39) 
(4.40) 



implies that t° 



an" 



(3° = a 2 U - P^i/U = 0. This with R = y/P - a^ 2 x /aU. 
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As expected, two of the eigenvalues are zero due to 
the fact that two rows of the matrix are zero. Of the re- 
maining eigenvalues, X^ 3 \ \^ and A^ are clearly real, 
but A± involve square roots through R so they require 
a more careful analysis. However, by inspecting the ex- 
pressions for A± it is easy to see that, though written in 
a different way, they are in fact identical to the eigen- 
values found in Ref. [1] for the case of the standard BM 
slicing condition. Since in that reference it was shown 
that those eigenvalues are always complex for f(a) < 0, 
always real for < f(a) < 1 and can always be made real 
for f(a) > 1 by an adequate orientation of the coordi- 
nate system, we conclude that our system of equations of 
equations is hyperbolic, at least weakly, as long as / > 0. 

The fact that two of the eigenvalues turn out to be 
identical to those found in the case of the standard BM 
slicing condition is surprising. One could have expected 
that, since the modified BM condition depends on the 
choice of shift, the shift vector should have affected these 
eigenvalues. The fact that it leaves these eigenvalues un- 
changed can be more easily understood in the 3+1 coor- 
dinate frame. We carry out this analysis in the Appendix. 

Having found that our equations are at least weakly 
hyperbolic, we now want to show that they are in fact 
also strongly hyperbolic. For this we must see if we have a 
complete set of eigenvectors. The eigenvectors associated 
with the matrix 1VP are: 



e± 



e (3) = 



[A ± ,l,0,0,(* x B ± ±ii) A±, 
^ y B±A±^ z B±A±] , 
1 (3 X \ n 



(4.41) 



1 



,0,0, 



H 1 



e (4) - [0,0,0,0,-5^,1,0], 
ii y/ x 

e(5) = [°' < ' 'nS%'°' 1] ' 

e< 6 > = [0,Z\1,0,Z 2 ,Z 3 ,Z 4 ] , 



; f , x \>,, U . \ \ . >!'• II , (4.42) 
\I/„\I>. 

(4.43) 
(4.44) 



,(7) 



[0,Z 5 ,0,l,Z b ,Z 7 ,Z«] 



(4.45) 
(4.46) 



where 



B± = 2(^A ± -1] , 



(4.47) 



a 2 U 2 



A± = 



(lP-*2)(P-a*|) 



(a- 2 - aV 2 x )Q* x 



(P- 2atf x W + (P v % + ^} , (4-48) 



W = 2n (/T - Q^ x ) . 



(4.49) 



The functions Z l are somewhat long expressions whose 
explicit form is not needed in what follows, so we will not 
write them here. 

It is well known that those eigenvectors corresponding 
to different eigenvalues are linearly independent. The 



eigenvectors associated with the degenerate eigenvalues 
A< 4 \A( 5 ) and A( 6 ),A( 7 \ can be seen to be independent 
from one another by inspection. Therefore, the eigen- 
vectors given above form a complete set. This, together 
with the fact that the eigenvalues are real for / > 
allows us to conclude that the system of evolution equa- 
tions (4.19,4.23,4.24) is strongly hyperbolic if / > 0. 



V. SINGULARITY AVOIDANCE 

In reference [1] it was shown that the original BM 
slicing condition can avoid so-called "focusing singulari- 
ties" [19] depending on the form that the function f(a) 
takes in the limit of small a. In particular, it was shown 
that if f(a) behaves as / = Aa n for small a and the 
normal volume elements vanish in terms of proper time 
r as 7^ 2 <~ (t s — r) m , one can have three different types 
of behavior: 

1. For n < the lapse vanishes before the normal 
spatial volume elements do, which corresponds to 
strong singularity avoidance 

2. For n — and mA > 1 the lapse vanishes with 
the normal spatial volume elements and the singu- 
larity is reached after an infinite coordinate time, 
corresponding to marginal singularity avoidance 

3. For both n > 0, and n — with mA < 1, the 
lapse vanishes with the normal volume elements but 
the singularity is still reached in a finite coordinate 
time, so there is no singularity avoidance. 

The results summarized above, however, depended 
crucially on the fact that the original BM slicing con- 
dition relates the evolution of the lapse to the evolution 
of the normal volume elements. The modified version of 
the BM condition, on the other hand, relates the evolu- 
tion of the lapse to the evolution of the coordinate volume 
elements. We must therefore see how this affects the con- 
clusions about singularity avoidance. The first thing to 
notice is that all the analysis of reference [1] will follow 
exactly in the same way for the modified BM condition if 
we replace normal volume elements with coordinate vol- 
ume elements. This means that the modified BM condi- 
tion will avoid singularities where the coordinate volume 
elements vanish (i.e. coordinate focusing singularities) 
under the same conditions as before. 

There is one very important difference between "nor- 
mal" focusing singularities and "coordinate" focusing sin- 
gularities. When the normal volume elements vanish, the 
normal direction to the hypersurfaces becomes ill-defined 
and the hypersurfaces become non-smooth. When the 
coordinate volume elements vanish, on the other hand, 
it is only the time lines that cross. This means that one 
could in principle develop a coordinate focusing singu- 
larity on a perfectly smooth hypersurface, or worse still, 
one could develop a normal focusing singularity for which 
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the time lines do not cross and the coordinate volume 
elements remain non-zero. The second case would be ex- 
tremely problematic as our hypersurfaces would become 
non-smooth but the lapse would not collapse in response 
to this. 

To see under what conditions we can have one type 
of focusing singularity and not the other we must look 
at the evolution equations for the normal and coordinate 
volume elements: 

dtln^J 2 = ~aK, (5.1) 
5 t ln 7 1 / 2 = -(aK-Dip 1 ) . (5.2) 

From the first of these equations it is clear that for a nor- 

l li 

mal focusing singularity to develop (-fj — + 0) we must 
have K — > oo. From the second equation we then see 
that the only way in which we can have a normal focus- 
ing singularity develop while at the same time keeping 
a non-zero coordinate volume element is for Di(j % to di- 
verge with K while keeping their difference finite. 

We then conclude that if the divergence of the shift 
remains finite, both types of focusing singularities will 
happen at the same time. This means that if the shift 
vector remains regular, the modified BM slicing condition 
will avoid singularities exactly in the same way in which 
the original condition did. 



VI. GAUGE SHOCKS 

In Section IV we have shown that the system of equa- 
tions (4.35) for the variables (4.34) is strongly hyperbolic. 
We can then define a complete set of "eigenfields" Wi in 
the following way: 



u = Ruj 



UJ = R~ 1 M 



(6.1) 



where R is the matrix of column eigenvectors ej. 

We say that the eigenficld uii is "linearly degenerate" 
if its corresponding eigenvalue Aj is independent of the 
eigenfield, that is 

IT = Y, 7T- 7T- = VuK ■ et = . 6.2 

3 = 1 

Linear degeneracy guarantees that the corresponding 
eigenfield will not develop shocks. 

Using the eigenvalues and eigenvectors found in 
Sec. IV, the conditions for linear degeneracy become 



C± := A± 



d\+ dX A 







n (3) d\W (3) d\W (3) d\& 
63 := " -- " 



dU 

(3) 9A(3) 

"T e 6 



d(3 x 



dp ' 7 d(3 z 
(4) dX^ d\W 



(6.3) 

(6.4) 
(6.5) 



C 5 
C 6 



J5) 



C 7 



df3 x 
= . 



d/3 z 



= 



(6.6) 
(6.7) 



A straightforward calculation shows that equations 
(6.4)-(6.6) are satisfied identically. On the other hand, 
equation (6.3) is precisely the same equation found for 
the standard BM slicing condition in reference [5] , where 
it was shown that it leads to the following condition on 
the function f(a) 



j 2 







(6.8) 



This means that the modified BM slicing condition is 
linearly degenerate under the same circumstances as the 
standard BM condition, and therefore it will avoid gauge 
shocks exactly in the same cases. 



VII. CONCLUSION 

We have studied a modified version of the BM slicing 
condition that has to important features: 1) it guarantees 
that if a spacetime is static or stationary, and one starts 
the evolution in a coordinate system in which the metric 
coefficients are already time independent, then they will 
remain time independent during the subsequent evolu- 
tion, and 2) the modified condition is naturally adapted 
to the use of a densitized lapse as a fundamental variable. 

By analyzing this modified BM condition written in 
covariant form on an arbitrary background spacetime, 
we have also shown that it is strongly hyperbolic for 
f(a) > 0, just as the original BM condition was. More- 
over, we have found that the characteristic speeds of the 
original BM condition are not modified. Finally, we have 
shown that as long as the shift vector remains regular, 
the modified BM condition avoids both focusing singu- 
larities and gauge shocks under the same conditions as 
the original BM condition. 

Because of these results we believe that the modified 
BM condition might be just as useful as the original BM 
condition for evolving strongly gravitating systems, while 
at the same time having the extra benefits described 
above. We plan to carry out numerical experiments to 
test this, but since the modified BM condition leads to a 
different slicing of spacetime for different choices of shift, 
such experiments will require first that one studies dif- 
ferent shift conditions. We are currently working on this 
issue, and will report our findings in a future work. 



Acknowledgments 

We thank Olivier Sarbach for useful discussions. This 
work was supported in part by CONACyT through the 
repatriation program and grants 149945, 32551-E and 
J32754-E, by D G APA-UN AM through grants IN112401 
and IN122002, and by DGEP-UNAM through a comple- 
mentary grant. 



10 



APPENDIX 

While studying the hyperbolicity of the modified BM 
slicing condition in Section IV we found, surprisingly, 
that two of the eigenvalues associated with this modified 
condition are identical with the eigenvalues associated 
with the standard BM condition found in Rcf. [1]. The 
analysis of Section IV was done in a covariant way, which 
required the introduction of an evolution equation for the 
shift that stated the fact that the shift has to be paral- 
lel to the hypersurfaces. The introduction of the shift 
equation made the analysis considerably more compli- 
cated and makes it difficult to see why the eigenvalues 
should remain equal. Here we will do a simple analy- 
sis in the 3+1 coordinate frame for the particular case 
of one spatial dimension to try to understand why the 
characteristic speeds are not modified. 

Let us then consider first the standard BM slicing con- 
dition in one spatial dimension. The evolution equation 
for the lapse has the form 



d t a = [3d x a - a 2 fK . 

If we now introduce the first order quantity A := 
we can rewrite this equation as the system 



(A.l) 



d t a 
d t A 



PA - a 2 fK , 
d x (pA - a 2 fK) 



(A.2) 
(A.3) 



with K the trace of the extrinsic curvature. On the other 
hand, from the ADM equations we find the following evo- 
lution equation for K 



d t K ~ d x <fiK - Ah) 



(A.4) 



where the symbol ~ denotes equal up to principal part, 
and where 7 is the one-dimensional coordinate volume 
element. Notice that we are considering the shift (3 to be 
a prescribed function of space and time. 

The Jacobian matrix associated with the evolution 
equations for (A, K) is then 



M 



-13 
l/ 7 



7 



(A.5) 



Notice that the overall sign chosen here is the one ob- 
tained by moving the spatial derivatives to the right hand 
side of the equations, as this is the sign we need if we want 
to associate the eigenvalues with characteristic speeds. 
The associated eigenvalues are 



A± = -/3±a(// 7 ) 1/2 
with corresponding eigenvectors 



v± 



±a{.f 1 ) 1/2 



(A.6) 



(A.7) 



We then see that evolution equations for the pair (A, K) 
form a strongly hyperbolic system as long as / > 0. 

Let us now consider the modified BM slicing condition. 
The evolution equation for the lapse now becomes 



3 t a = -af [aK - 8 x p - (/3/2 7 ) 8 xl ] 



(A.8) 



Notice how this now includes a spatial derivative of the 
coordinate volume element, which can not be considered 
as a source since by construction the time derivative of 
the lapse is proportional to the time derivative of the 
volume element. If we now define G := d x -f, we can 
rewrite the evolution equations for a and 7 as the system 



d t a 



-a 2 fK + 



afP 
2 7 



G + afd x p , 



d a = ~2ajK + 0G + 2jd x p 
d t A ~ 8 X (-a 2 fK 
8 t G ~ d x (-2ajK + pG) . 



(A.9) 
(A.10) 
(A.11) 
(A.12) 



We now see that the Jacobian matrix associated with the 
system (A, G, K) takes the form 



M 




-afp/2 7 

-P 




a 2 f\ 
2a-/ 



(A.13) 



d x a, The eigenvalues of this matrix are easily found to be 



A = 0, 

A± = -p±a(f/ 1 ) 1/2 



with associated eigenvectors 



''0 



(P 2 ll 2a 1 ,p) , 
= (±a(/7) 1/2 ,±2 7 3/ 7/ 1/2 ,l) • 



(A.14) 
(A.15) 



(A.16) 
(A.17) 



We can again see that two of the eigenvalues correspond- 
ing to the modified BM slicing condition coincide with 
the eigenvalues of the standard BM condition. However, 
the reason for this is now much easier to see. Notice 
that the Jacobian matrix (A.13) has two rows that are 
multiples of each other (row two can be obtained from 
row one by multiplying with 2j/af). This means that 
we can define a new variable £ := G — (2j/af)A that 
will evolve only through lower order terms. We can now 
make a change of variables from (A, G, K) to (A, E, K), 
by replacing G with S + (2 r y/af)A. Since E evolves 
only through lower order terms, its derivatives can be 
treated as source terms in the hyperbolicity analysis. It 
is then easy to see that the Jacobian matrix for the re- 
duced system (A, K) is identical to the Jacobian matrix 
associated with the standard BM slicing condition given 
in Eq. (A.5). 

The important observation is that even though the 
modified BM slicing condition replaces the Lie deriva- 
tive of a with respect to the shift with the divergence 
of the shift, since this divergence includes the Lie deriva- 
tive of 7, and the time derivative of a and 7 are multiples 

of each other, we recover in the end precisely the same 

1/2 

characteristic speeds —f3 ± a (f/j) . 
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